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The excitation spectrum and the band structure of a Bose-Einstein condensate in a periodic 
potential are investigated. Analyses within full 3D systems, finite ID systems, and ideal periodic 
ID systems are compared. We find two branches of excitations in the spectra of the finite fD 
model. The band structures for the first and (part of) the second band are compared between a 
finite ID and the fully periodic ID systems, utilizing a new definition of a effective wavenumber 
and a phase-slip number. The upper and lower edges of the first gap coincide well between the 
two cases. The remaining difference is explained by the existence of the two branches due to 
the finite-size effect. 
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1. Introduction 

Bose-Einstein condensate (BEC) of the atomic gases 
trapped in the periodic potential of an optical lattice 
are currently intensively investigated.^""^ These systems 
have been researched also as the candidates of quantum 
computing.'^' ^ The condensate is usually confined with a 
potential generated by a magnetic field, and also mod- 
ulated by the periodic potential formed by the counter- 
propagating laser beams; for a schematic illustration of 
the experimental arrangement, see Fig. 1. 

The system may be treated as periodic^' ^ ignoring 
the confinement. Within this periodic framework, the 
wavenumber of the excitations becomes periodic and the 
whole spectrum splits into band structures which are 
similar to those within solid-state physics. Conversely, 
in the limit of a weak periodic potential, the conden- 
sate reduces into a harmonically trapped, cigar-shaped 
condensate. The classifications of the collective excita- 
tions are closely related to the symmetry of the confining 
trap. Angular momenta and other indices in spherically 
symmetric systems have been used.^ Within cigar-shaped 
systems, the wavenumber along the z-axis* is also a use- 
ful quantum number to characterize an excitation. 

In this paper, we calculate the full excitation spectra 
for Bose-Einstein condensates of atomic gases trapped 
in a periodic optical potential and a harmonic confine- 
ment potential. We aim to classify and understand the 
excitation spectra properly, using the similarity with the 
purely periodic systems and the cigar-shaped systems. 
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Fig. 1. BEC of Rb atoms is confined along the z and r directions 
by harmonic traps with frequencies and Ur respectively. The 
condensate consists of 5 X 10'' atoms of *^Rb. Counterpropa- 
gating laser beams of wavelength A form a potential of period 
d= A/2. 



2. Condensate in ID Optical Lattices 

We study the excitation spectra of Bose-Einstein con- 
densates (BEC) confined in a composite potential 
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that consists of an optical lattice in the principal direc- 
tion of consideration and harmonic confinement poten- 
tials along the lateral dimensions. Here s is an inten- 
sity parameter of the periodic optical potential, i?R = 
^ (^) is a recoil energy, to is the mass of a ^''Rb 
atom, Vr and are the frequencies of the harmonic po- 
tentials, and d is the period of the potential. We employ 
d = 0.395 /im which is one-half of the typicaP wavelength 
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A for the laser beams, Vr ~ 100 s~^ and I'z = 30 s~^ 
which are also within the typical ranges in experiments. 
A condensate with 3D geometry, ID geometry and a pe- 
riodic ID situation arc simulated and compared in order 
to determine the effects of the confinement potentials. 

The BEC in an optical lattice with the lower value of 
the intensity s is well described by the Gross-Pitaevskii 
equation^ 

(-CV2 + y(r, z) + g\^{r, z)\^) 0(r, z) = ^l<f>{r, z) (4) 

where (j){r, z) is the wavefunction of the condensate, /i is 
the chemical potential, g = ^'^^ ° is an interaction pa- 
rameter and C = fi? /{2m) is constant. We treat a con- 
densate with 5 X 10^ atoms of *^Rb. The scattering length 
a = 5.4 nm is employed for the ^^Rb atoms. The period 
of the optical lattice is = 0.395 /xm and one lattice site 
is occupied by at most 1300 atoms. The condensate oc- 
cupies about 80 sites. Figure 2 shows the density of 
a typical condensate along the ^;-axis. 
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Fig. 2. Solid curve: the condensate density \ij>{r,z)\'^ at r = 
with the lattice intensity s = 5. The range z > is plotted. The 
dotted line is the density profile in the absence of the optical 
lattice (s = 0). 



The condensate supports excitation spectra within the 
Bogoliubov framework. The Bogoliubov equations'' 

{-CV^ + V{r,z) + 2g\,i)\^ - H)uq-g(fvg = SqUq, 

(5) 

{-CV^ + V{r,z) + 2g\<i>\^ - n)vq-g<i>*^Uq = -SgVg, 

(6) 

yield the excitation energies Sq and the corresponding 
wavefunctions, Ug and Vg. The condensate and the exci- 
tations may have finite angular momenta which are es- 
sential to treat the systems having vortices.^ We only 
consider excitations having zero angular momentum and 
concentrate on the wavenumbers along the principal z- 
axis in this article. 

3. Wavenumbers 

In the absence of an optical lattice potential, the con- 
densate has a cigar-like shape and the wavenumber of 
an excitation may be defined through differentiating the 
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Fig. 3. Excitation spectra within 3D analysis (crosses), finite-size 
ID model (bullets) and periodic ID analysis (solid line), (a) The 
wavenumber k defined by Eq. (8); and (b) the cfTcctive wavenum- 
ber fcgfi as given in Eq. (9). Both ID spectra (bullets) feature two 
distinct branches, that are named branch A and B. The patterns 
of the 3D spectra significantly differ between the right- and left- 
hand sides for branch A. The range of feefi does not depend on s. 
The vertical lines in (a) and (b) indicate fe" and k^g respectively. 



wavefunction. In contrast, we may consider an ideal peri- 
odic ID system to treat only the effects of the ID optical 
lattice potential and ignore the confinement potential. 
Within the periodic ID system,^ the wavefunction may 
be expressed in the form 



= eiip{iKz)u^{z), 



(7) 



where k is a quasi- wavenumber. 

In order to consider these two aspects of this system, 
we introduce the two definitions for wavenumbers 
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and 



fceff(<i) 



, .12 
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where Uj is the vahie of the wavefunction u(r, z) at the 
center of the j-th site. The definition of kes arises from 
the discretization of the second derivative in eq. (8). 
Hence the two wavcnumbcrs k and k^ff coincide in the 
limit as d ^ 0. The keS is compatible with the quasi- 
wavenumber k within a periodic system for \k\ -C 1/d. 
The corresponding wavenumbers of the condensate are 
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l^eS — 
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Figures 3(a) and 3(b) represent the excitation spectra 
using k and keS, respectively. Note, in particular, that the 
wavenumbers k are strongly shifted and modified with 
respect to those (|k| < ir/d) of an infinite periodic ID 
system.^ The wavenumbers k are appreciably high com- 
pared to w/d. This is because the intensity of the wave- 
function u oscillates with the period d along the length 
z. The effective wavenumber k^s defined in eq. (9) is not 
explicitly modified by this oscillation. Its value does not 
exceed ir/d. The lowest fccff value is close to 0, which is 
the minimum |k| in the purely periodic system. 

The alignments of the excitations in the plots become 
disordered [Fig. 3(b)] and dense [Fig. 3(b)] above certain 
threshold values of k and keS. The spectra (the bullets 
in Figs. 3(a) and 3(b)) of the corresponding ID model 
systems make the change in the alignments plausible. 
The ID model system is calculated by treating only the 
r = case in eqs. (1-6). The chemical potential /it is the 
same as that of a corresponding 3D system, such that 
the particle densities are close to those of a 3D system 
at r = 0. 

The ID system features two branches in the spectra, 
which are named A and B. The spectra in the 3D anal- 
ysis differ considerably between the left- and the right- 
hand sides of the branch A in Fig. 3(b). The branch A 
explains the dense population of the excitations around 
k = O.dir/d (hidden behind the bullets in Fig. 3(a)) and 
fceff = 0.57r/(i in the spectra of the 3D system. The 3D 
systems have modes corresponding to the branch A with 
various radial wavenumbers in the narrow range of the 
wavenumber fccff of principal direction. 

The vertical (i.e., having definite k and fcefj values) 
branch A does not exist within the purely periodic model. 
For comparison, we consider spatially periodic systems. "'^ 
The condensate wavefunction has the same form as in 



eq. (7) and Uk{z) obeys the wave equation 
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+ (fl'l^i«(^)P + Vopt) u^{z) = ij,{k)u^{z) (12) 

where — ^ < k < ^ . The potential energy V^pt is the last 
term in eq. (1). The period of the system is d and the av- 
erage density of atoms is 10^° m~^, which is comparable 
to the densities within the previous 3D and ID analyses. 

The excitation energies are written using the energy of 
a condensate wavefunction 

= I -Cut { — Ik I Uk{z) 
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(13) 



The solid line in Fig. 3(b) represents the spectra = 
{E^ — Eq)/N, which are comparable to the previous spec- 
tra that utilize fceft- 

4. Origin of the Branches 

Figure 3(b) shows that branch A does not exist in the 
spectra of the purely periodic system. In order to explain 
the origin of branch A, we investigate the profiles of the 
wavefunctions within the ID system. Figure 4(a) shows 
the spectral evolution within the ID analyses. Excita- 
tions in the branch A are localized outside of the con- 
densate \z\ > Ztf- Here ztf = 14.6 /xm is one-half the 
condensate length within the Thomas-Fermi approxima- 
tion in the absence of the periodic potential. The Bogoli- 
ubov eqs. (5) and (6) reduce to a Schrodinger equation 
outside the condensate (|^| — > 0). The graph ofVz{zpeak) 
in Fig. 4(b) shows that the increase of the excitation en- 
ergies in branch A simply depends on the z-confinement 
potential T4(-2peak), where Zp^^ak indicates the peak in- 
tensity of an excitation wavefunction tiq{z). The regular 
spacing of the excitation energies e may be understood 
in terms of the ^-confinement trap and the period d of 
the optical lattice. The energy spacing Ae reflects the 
z-confinement potential since 



Aec^V{z + d)-V{d). 



(14) 



Excitations within the 3D framework in Fig. 4(b) have 
the same features. Branch A in the ID systems corre- 
sponds to excitations localized outside the condensate in 
3D systems. 

The origin of the branch B is clearly seen through 

defining a number / of phase slips. Unlike for the purely 
periodic system, the differences in the complex phases of 
a wavefunction Uq{z) between the neighboring sites are 
always either or tt. Therefore, we may deflne a phase- 
slip number 
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Fig. 4. (a) Intensities of the excitations |?t(2)p + |t)(z)p within 
the ID framework. There is a sijatial separation between the 
excitations of the two branches A and B in Figs. 3(a) and 3(b). 
(b) Excitation energy vs potential energy Vz{z) = ^{2ttUz)'^ z'^ 
at the peak values of the amplitude of wavefunctions \u\^. The 
relation is linear beyond (zpeak > ztf) the condensate, both in 
the 3D and within the ID framework. Here ztf = 14.65 jum is a 
half the Thomas-Fermi length of the condensate for s = 0. The 
vertical line indicates Vz{ztf)- The slope of the solid oblique 
line is 1. (c) The number of the phase slips /. Labels A and B 
indicate the two branches. The smooth increase in / breaks down 
at the branch separation. The intensity s = 5 was used in these 
computations. 



for each of the excitations. Figure 4(c) shows this number 
/ as a function of the effective wavenumber, k^s- The 
phase-shp number / increases as kes increases from zero 
towards the branch separation. In other words, the peak 
of the phase-shp number separates the branch B from 
the rest of the excitations. 

5. Band Gap 

The ID systems feature gaps in the spectra if we ignore 
the branch A as shown in Figure 5(a). The range of the 
gap is close to the first gap in the corresponding periodic 
system. The gaps between the first and the second bands 
are plotted in Fig. 5(b) as functions of s for the range 
1 < s < 10 of intensities for the optical lattice. The upper 
limit of the gap, i.e., the upper solid line and the upper 
scries of bullets in Fig. 5(b) coincides between the two 
systems. There is difference of about 0.17£^h between the 
two lower limits, which are the lower solid line and the 
lower series of bullets in Fig. 5(b). 

The excitation energy of a wavefunction which has the 
maximum number / of phase slips indicates the separa- 
tion of the branches. The open circles in Fig. 5(b) show 
the excitation energies at the separation. These energies 
coincide better with the lower limit of the band in the 
periodic system. The difference between the two lower 
limits emerges from the existence of the branch B. There- 
fore, the two branches make a significant difference in the 
spectra between a finite system (ID) and an ideal peri- 
odic system. Otherwise, the two spectra feature quite 
similar first bands and first band gaps. 

6. Summary 

We have investigated the excitation spectra and the 
band structures of a Bose-Einstein condensate in periodic 
potentials. Comparisons between calculations performed 
within 3D systems, finite ID systems, and ideal periodic 
ID systems have been made. 

Two branches are found in the spectra of the ID model 
system. The branch A, which makes a major difference 
between the confined ID system and a periodic ID sys- 
tem, arises from the excitations that prevail outside the 
condensate. Excitations in the other branch B possess 
lower phase-slip numbers than the excitations at the 
branch separation, despite the higher wavenumbers k and 
fceff- The structures corresponding to the branches are 
found in the spectra of the 3D system in Fig. 4(b). 

The band structure up to the beginning and along the 
2nd band are compared between the finite ID system and 
the ideal periodic ID system, utilizing new definitions of 
a effective wavenumber and a phase-slip number. The up- 
per and lower edges of the first gap coincide well between 
the two systems. The remaining differences are explained 
by the existence of branch B in the excitation spectra of 
the ID system. 
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Fig. 5. (a) Excitation spectra of the ID system and tlie periodic 
system for the optical lattice intensities s = 1 and 2. Both spectra 
display a band gap indicated by the arrows, (b) The gap between 
the first and second bands is plotted as a function of s. The solid 
lines are the upper and lower edges of the gap within the periodic 
system. The bullets are those for ID systems. The open circles 
show the excitation energies of wavefunctions which have the 
maximum number / of phase-slips for each s. 



